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Origin of black string instability 
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It is argued that many nonextremal black branes exhibit a classical Gregory-Laflamme (GL) 
instability. Why does the universal instability exist? To find an answer to this question and explore 
other possible instabilities, we study stability of black strings for all possible types of gravitational 
perturbation. The perturbations are classified into tensor-, vector-, and scalar-types, according to 
their behavior on the spherical section of the background metric. The vector and scalar perturbations 
have exceptional multipole moments, and we have paid particular attention to them. It is shown 
that for each type of perturbations there is no normalizable negative (unstable) modes, apart from 
the exceptional mode known as s-wave perturbation which is exactly the GL mode. We discuss the 
^vQ , origin of instability and comment on the implication for the correlated-stability conjecture. 
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I. INTRODUCTION 



Stability of a given spacetime is a crucial issue from many standpoints. In general relativity, a stable spacetime 

■ will be realized by a dynamical evolution starting from a generic set of initial data on a Cauchy surface. However 
stability in general relativity is frequently subtle issue, and because of that it becomes important and interesting in its 

£NJ ■ own right. From a string theory perspective, it is interesting to know what spacctimes are appropriate backgrounds 
for studying string propagation and its dynamics. Besides, information of gravitational dynamics and properties are 
useful to understand Yang-Mills theory by means of gauge/gravity dualities, and vice versa In 
this respect, instability on the gravitational side is an indicator of interesting gauge theory dynamics, such as phase 
transition and so on 

The fundamental generic instability is Gregory-Laflamme (GL) instability |lC( , which is accompanied by a uniformly 
Cl^, smeared horizon. The fundamental phenomena is however one of long-standing puzzles in gravity. For example, (i) 

■ what is the necessary and sufficient condition for the onset of a dynamical instability of a horizon ? (ii) Why is a uni- 
form horizon unstable? The first question was addressed by a so-called correlated-stability conjecture (CSC) [TlllT^|. 
Namely, the onset of the dynamical instability of black brane will be the same as the onset of (local) thermodynamic 
instability. The second question is more fundamental and naive. The origin of the instability might have deep con- 
nection with quantum aspect of gravity, since the onset of instability is predictable by black hole thermodynamics 
due to CSC. Here we would like to pursue the question from classical aspect of gravity. (See 01 f° r fluid analogy of 
GL phenomena.) 

First of all, we do not know full dynamics of unstable black objects in higher dimensions |l4j| . In particular, as far 
as the present author knows, a complete analysis of (in)stability has not been carried out. (A numerical investigation 
for the 5-dimensional black string in the braneworld model with AdS bulk was performed in (l 5| - ^ In fact, even for 
the higher dimensional Schwarzschild black holes (BHs), its dynamical stability was established in recent years by 
Kodama and Ishibashi (KI) 0, 0, . The instability found by Gregory and Laflamme is the s-wave mode, and 
the perturbation is "minimum" deformation of horizon. For perturbations with higher multipole moments, similar 
instability might persist. An interesting point is that existence of instability implies existence of a critical static mode 
and the mode could be continued to a state with nonperturbatively deformed horizon [T^. l2fl l2ll 12^ . l23l | , so that 
any extra instability implies extra static sequence of solutions. Besides, they will have physical meaning in Euclidean 
space |2J, I23, |2^, I23 . 1 In addition to the stability issue of uncharged black branes, complete stability of BPS state 
with respect to all possible types of perturbations, which should include breaking of supersymmctry, remains an open 
question, although there are several evidence for it 

In order to promote greater understanding of the nature of black string/brane, it is inevitable to investigate the 
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1 The stability argument in Ref. L'4|' is sometimes applied to black strings, using Wick rotation. The proof of stability for higher 
multipole moments assumes that all eigenvalues are real under the periodic Euclidean time. In general, this assumption for eigenvalues 
and boundary conditions is crucial for stability argument, and we should not naively apply the argument to discuss the stability of black 
strings. 
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stability with respect to all the types of perturbations. Following to the general gauge- invariant formalism for higher 
dimensional maximally symmetric BHs by KI, we develop a general perturbation theory of black string and tackle 
the stability problem. (See also [3ll I32I ] for the basic work related to the gauge-invariant formalism of maximally 
symmetric spacetimes.) In this approach, the perturbation variables are classified into three types, those of tensor, 
vector and scalar modes, according to the type of harmonic tensor used to expand the perturbation variables. Contrary 
to the perturbations for the maximally symmetric BHs, vector and scalar type perturbations will not have simple 
master variables due to extra physical degrees of freedom. We study stability of these perturbation variables. 

The paper is organized as follows. In the next section we first classify perturbations into tensor, vector, and scalar- 
types with respect to the maximally symmetric 71-dimensional spacetimes. Then for each type of perturbations, we 
express the Einstein equations in terms of them. In Sec. IIIII we study stability of tensor and vector perturbations 
and found that there is no instability in these perturbations. In Sec. IIVI the stability analysis of scalar perturbation 
will be carried out. The unstable GL mode is an exceptional mode in the present perturbation scheme and we discuss 
that there is no other unstable mode in the black string perturbations. The origin of such an exceptional mode will be 
clarified in comparison with the perturbations for the maximally symmetric BHs. Section is devoted to summary 
and discussion. Throughout this paper we follow the notation in Refs. p^Hr^ ll^. 



II. GENERAL PERTURBATION THEORY 



As our background spacetimes, we consider the D = (n + 3)-dimcnsional metric of the form 

ds% = g A Bdx A dx B = g ab dy a dy b + r 2 da 2 n + dz 2 , (1) 

where g a b is the Lorentzian metric of the two-dimensional orbit spacetimc, and da^ — r )ij{y)dy l dy^ is the metric of the 
n-dimcnsional maximally symmetric space K n with sectional curvature K — 0, ±1. Throughout this paper, we use 
the notation a, b = 0, 1, i, j = 2, • ■ ■ , n + 1 and a, j3 ■ ■ ■ = 0, ■ • ■ , n + 1. The covariant derivative with respect to the 
metric g a b and 7^ arc defined as D a and Di, respectively. In the followings, we develop a general perturbation scheme 
for black objects with co-dimension one. The perturbation will be specialized to the black string perturbations in the 
next section (See footnote |21 for more general perturbations.) 

Most general metric perturbations SgAB for this background spacetimes are 

ds 2 = {gAB + 8gAB)dx A dx B . (2) 

Utilizing gauge degrees of freedom, xa — * x' A = xa + £a(%p), we can eliminate perturbations in z direction at any 
times, taking a Gaussian normal coordinates: 

ds 2 = (g a p + 5g a p)dx a dx 13 + dz 2 . (3) 

This gauge fixing is however not complete. There are two types of residual gauge degrees of freedom. The corre- 
sponding infinitesimal coordinate transformations are 

£ z = P(x a ), & = -zd a P(xp), (4) 

and 

C = 0, £ a = Q a (xp). (5) 

The first one corresponds to shifting z — const, surface. The second is the gauge transformation transverse to a 
z = const, surface, and hereafter we call this "transverse" gauge degrees of freedom. 

Because the background spacetimes are translationally invariant along z direction, we can take arbitrary hypersur- 
face of 

z = const. = z c (6) 

to study the perturbations without loss of generality. This approach is a sort of an effective theory approach. In this 
approach, the residual gauge P is fixed once we take a z = const, surface, on which we will study perturbations. 

At this point, if we consider only homogeneous perturbations along the z direction, the general perturbations J3J is 
the same as the gravitational perturbations of maximally symmetric black holes in higher dimensions studied by Ko- 
dama and Ishibashi. Their perturbation theory is most generic and based on gauge-invariant scheme, yielding master 
variables for each type of perturbations. Following to their perturbation theory, below, we develop transversely gauge 
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invariant perturbation theory, so that perturbation variables independent of the residual gauge JJJJ are introduced. 
A point is that the general perturbation provides transparent perturbation scheme, which can be directly compared 
with the perturbations of maximally symmetric black holes. 2 

A. Tensor-type perturbations 

We begin by considering tensor perturbations, which are given by 

Sg a b = 0, Sg ai = 0, 8g^ ; = 2r 2 if T T.y , Sg az = 0, (7) 

where Ht is a function of {£, r, z}, and the harmonics tensors Tij are defined as solutions to the eigenvalue problem 
on the n-sphere; 

(A„ + fc 2 )Ty = 0, T*< = 0, = 0. (8) 

Here, Dj is the covariant derivative with respect to the metric 7^ and A„ = r ) lJ DiDj. In these equations we have 
omitted the indexes labeling the harmonics and the summation over them. For K = 1, the positive eigenvalue k\ for 
a discrete set, k\ = 1(1 + n — 1) — 2, / = 1, 2, • ■ ■ . 

The tensor perturbations are essentially transversely gauge invariant. Following [Til IT7L ll8L l3lj . we introduce a 
new variable <& = r"/ 2 iJ-r- The master equation follows from the vacuum Einstein equations 1B13() . 



□ -i 



fcy + 2?iA H rUr H (Dr) 2 



* + *,» = 0, (9) 



where □ = D a D a denotes D' Alembertian operator in the two-dimensional orbit space. We remind that there is no 
tensor-type harmonics on a 2-spherc, so that the tensor perturbations only exist for n > 3. 

B. Vector-type perturbations 

Perturbations of the vector type can be expanded in terms of vector-type harmonic tensors Vj satisfying 

(A„ + A$)Vi = 0, l) s " 0. (10) 

As in the case of tensor-type harmonics, the eigenvalues ky are positive definite. For K = 1 the eigenvalues form a 
discrete spectrum given by 

kl = 1(1 + n- 1) - 1, l=l,2,---. (K=l) (11) 
In terms of vector harmonics, metric perturbations are expanded as 

Sg ab = 0, Sg al = rf a Vi, Sg l3 = 2r 2 H T Y lv Sg az = 0, (12) 

where and Vj satisfy 

Vij = -^-(DiVj + DjVi), V% = 0, D j W i = kv ~^~ 1)K V i . (13) 

ZKy zfcy 



2 We have decomposed the metric into 2 X n X 1 space with employing the gauge fixing. More general formulation will be possible by 
decomposing the metric into m X n space 13 ll , where m is m > 3 depending on the translationally invariant spatial dimensions of black 
bra ne. By employing such decomposition, we can use many covariant formulas for the higher dimensional maximally symmetric BHs in 
1 1(1 without significant changes, although such fully gauge invariant equations give more messy equations of motion for each variables. 
In this picture, it is easy to count a number of physical degrees of freedom for each type of perturbation. The physical degrees of freedom 
for tensor and vector F a (see Eg. 1141 ^. are 1 X [Tij] and (m — 1) X [Vj], respectively, taking into account the number of constraint 
equations for vector perturbation. Here [Tij] = (n + l)(n — 2)/2 and [Vi] = (n — 1) are the number of degrees of freedom for the respective 
harmonics. The physical degrees of freedom for scalar perturbation, F and F%, are (m 2 + m + 2)/2 — (m + 1) = m(m — l)/2, subtracting 
the number of constraint equations for scalar perturbation. The total gravitational degrees of freedom are (n + m)(m + n — 3)/2. 
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Note that also satisfy [A + k v — (n + l)K]Yij = 0. The special mode k v = (n — 1)K is known as the exceptional 
mode for the vector perturbations, since V,j vanishes for this mode. 
For k v 7^ (n — 1)K > 0, transversely gauge invariant quantity is 



F a (t,r,z) = fa + rDa 



Ht 
ky 



The vacuum Einstein equations, SRij = and 5R ai = 0, reduces to 



d 2 (Ht 
dz^ \~k~v~ J ' 
„ (r n+1 F C 1 )) - m v r n - 1 e ab F b = -r n+1 e ab f b 



D a {r n - 1 F a ) = -r n 
D 



where my = k v — (n — 1)K , and we have introduced 



r e ab D 



(14) 

(15) 
(16) 

(17) 



The Einstein equation SRi Z = gives a non-vanishing equation, but it is not an independent equation. Combining 
these two equation, we obtain an equation for F a , 



ad 7-> 


Y^7) 


+ r n+2 D a 









^D^F C ) 



m v r n - 1 F a = -r n+1 F a zz . 



(18) 



Therefore our stability problem is reduced to solve the equation of motion (EOM) for the vector F a . The vector <|14[) 
has been constructed to be invariant under the gauge transformation which is independent of z. Thus any solutions 
of the evolution equation (|18|) have physical meaning. 

Here we note that for the zero mode d z d z F a = the divergenceless condition (|15f) holds for the vector F a . From 
this condition, a master variable can be introduced, and the second equation Htj|) with employing the master variable 
reduces the Regge- Wheeler equation for n = 2, K = 1. By contrast with the zero mode, the KK modes have one 
extra physical degree of freedom. The two physical modes are governed by Eq. l|18f) . which will give coupled second 
order differential equations. 

The exceptional mode k v = (n — 1)K , corresponding to K = 1 and I = 1, receives special consideration. In this 
case the perturbations variable Ht does not exist because V,j vanishes, and correspondingly, Eq. ijT5|l does not exist. 
F a is not invariant under the transverse gauge transformation, and it has only one physical degree. Taking Ht = 
in (|14|) . the single physical mode which is invariant under the transverse gauge is given by (|17|) . For the zero mode, 
the equation for F^ is D a {r n+1 F^) = 0, and its solution is F^ 1 ' = const. /r^ n+1 \ This solution corresponds to 
adding a rotation to the background solution, although it is not a dynamical freedom. For the massive spectrum of 
this exceptional mode, the transversely gauge invariant equation is from (|16|) 



D' 



1 



j.n-\-2 ' 



D c (r n+1 F^ 



.1^(1) 



(K =!,£=!) 



(19) 



C. Scalar-type perturbations 

Scalar perturbations are given by 

Sdab = fab§, 5g al = rf a Si, Sg i:j = 2r 2 (H L jijS + Ht^), Sg az = 0, (20) 
where the scalar harmonics S, the associated scalar harmonic vector Sj, and the traceless tensor Sy are defined by 



1 



(A n + fc|)§ = 0, Si = -— A§, 



1 



- 1,2 



1 



J D iJ D i S+-7ii§, 



(21) 



k s og 

with the eigenvalues kg given by fc| = 1(1 + n — 1) for K = 1. By definition, and §y have the following property 

',.2 
<S 



D l S t = fc s S, S*i = 0, D'Sij = 



n k s 



-S; 



(22) 
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We introduce F(t, r, z) and F ab (t, r, z) defined by 



F = H L + -H T + -D a rX a 
n r 

F a b = fab + D a Xb + D X a , 



(23) 

where X a = (f a + j^D a HTj • By using these expansions, we have calculated Einstein equations for the scalar 
perturbations, which are summarized in the Appendix IU1 

Let us first consider the equations for the generic modes kg (kg — nK) ^ 0. The equations directly obtained from 
the Einstein equations contain such as f ab z and Hl, zz as we see in Eqs. I|C3|) and (|C4j) . Eliminating such terms by 
utilizing I|C1J| and (|C2|) . we obtain the following perturbation equations for F a b and F: 



D a D b F ab - UF C C - n —(D a F c c - 2D b F ab ) 



D a D b r 



+ (n 2 - 3?i-2) — 



D n rD h r 



F 



ah 



+ '-§■ F c c + 2DF - [Ak'g - 2(n - l)(n + 2)K + 4(n + l)(n - 2)(Drf\ % + 2(n + 1)F,„ = 0, 



(24) 



£> c r 
r 



□F afc + — {nD c F ab - AD [a F b)c ) - 2R^ T 



2R ( ?LF cd 



k 2 

-§F ab + 2(71 - 2) 



F c(a D b) D c r F c[a D b) rD c r 



' D a D b r D a rD b r\ 
-D (a rD b) F - 4(n - 2) ( — -g— I F 



+ ( ^cd +n(n- 1) 



D c rD d r , 



F c c - 2tiDF 



g a b 
n + 1 

2n(n+ 1) 



i^r/F^ - DF C C - -D c r(D c Ff - 2D d F cd ) 



D c rD c F + 2(n — 1)- 



nK 



-F 



+F ab 



0, 



(25) 



where (a 6) is a notation for the totally symmetric parts of tensors |3^|. For the zero mode, I|C1(I and l|C2f> work as 
"constraint" equations. In the present case, they constitute d 2 X a , which is given by IjCllJI . 
Additional EOMs are obtained from 5Ra z = 0. From HC6f) and l|C7|l. we get 



d 2 jD b Fab 



2n^F - D b D a X b - DX a + n— ( 2^X C 
r r V r 



-F a 



D a X c - D c X a ) + -fX c 



+klr 2 D a 



1 / F c 



2 D a r 



F zz + -^—D b (r n - 2 F b ) + - -"^ Z*l - (n - 2) 



= 

D a D b {r n - 2 F b ) 



(n-1) 



(Drf 



□r 
r 



2 r n-2 

k 2 s - (n - l)K 



□r , , (Dr) 2 
r r z 



D a r { F c 



nF 



f<: 



+ (n-2)F 



□ 



F c 



+ (n- 1)F 







(26) 



(27) 



where (ICllll is used to calculate X a , zz - Finally, 8R ZZ = gives Eq. |(jj£ 

af(/ c c + 2nH L ) = 



(28) 

Substituting ilC8|) and (|C9|) into this equation, we obtain an equation which does not contain z-derivatives, in contrast 
to the above four equations. These five equations are the basic equations for I > 2 modes. We will analyze these in 
the next section. 

For the exceptional mode k 2 s (k 2 g — nK) = 0, we need special consideration for the metric perturbations since some 
harmonic functions vanish in this case. For k 2 s = nK, which corresponds to I = 1, Eq. 1C2(I does not exist since 
vanishes and Ht is not defined. For k 2 s = 0, which corresponds to I = 0, both Eqs. 1C1|) and (|C2|) do not appear 
since Ht and f a do not exist. In the following we will consider these two exceptional modes separately. 



For £ = 1 (k s = nK), the metric perturbation Ht and hence Eq. 1C2(I does not exist. In this case, Eq. I|23|) is 
replaced by just setting Ht = 0. The transverse gauge transformation of F and F ab becomes 



SF 



h 2 f T 

%L + D a rD, 1 



6 



8F ab 



-D„ 



r 2 ( L 
k s \r 



D h 



r 2 ( L 
k s \r 



(29) 



and they are no longer transversely gauge invariant. We will use this gauge degree of freedom when we explicitly 
solve this mode. 

Equations for F and F ab are obtained as follows. From Eqs. (|C5|) . I|C6|) and (IC7II 



i) 2 



D c r 

2F - D C X C -(n+ 2) X c 

r 

F c c +(n- 2)D C X C + n 2 —X c 

r 



0, 



0, 



D b F ab + n—F ac - n 
r 



,D a rD c r 



D a r 



X r - n-^D°X 



Here X a is given by (|C1|1 : 



■UX a - D"D a X h +[—) X a - n —(D a X c + D c X a ) 
r J r 



X a , zz = -^D b {r n ~ 2 F h a ) + rD a H± ) 



0. 



(30a) 
(30b) 

(30c) 
(31) 



Utilizing this X a , Eqs. l|C4f) (or (|C8|> ) and (|C10() can be written in terms of F and F ab . These five equations are the 
basic equations for 1=1 mode. 
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For the s-wave (£ = 0) perturbation, H T and f a do not exist since S, ; and Sy cannot be defined for this mode. 
Hence F ab and F are given by F ab = f ab and F = Hl- The equations for these variables are given by (jC3f) and JT3J). 
[or equivalcntly, (|C8|) and HC10|) ]. Other complementary equations are from SR az = and SR ZZ = 0, i.e., Eqs. (IC5I) 
and P))l . 



d 2 z {F c c + 2nF) = 0, 

)2 



[D c F^ + n^F ca -2n^fF\ 



The variables F ab and F are gauge dependent, and their four components are reduced to two physical degrees of 
freedom by fixing transverse gauge (on z =const.). For example, the harmonic gauge condition, V a 1iab = 0, is a 
useful gauge fixing, which gives 



D c r D c r 

D F ac + n F ac - 2n F = 0. 

r r 



(32) 



We will use this gauge fixing later. 



III. STABILITY ANALYSIS: TENSOR AND VECTOR 

A. Background spacetimes and stability condition 

In this section, we discuss about stability of a higher dimensional black string. As a black string solution, we 
consider a following metric form: 



g ab dx a dx b = -f{r)dt 2 + 



dr 2 
W) 



(33) 



f(r) = K - 



n-l 



7 



Here g a b is the Lorentzian metric of the two-dimensional orbit spacetime, as mentioned in the previous section and the 
constant parameter r% defines the horizon radius. Hereafter, we only focus on the K = 1 case, because our interest is 
in the stability of the black string whose intersect is the higher dimcnsionally Schwarzschild black hole. 3 

If the equations of perturbations are reduced to a 2nd-order Schrodinger-type differential equation, the analysis of 
the stability can be carried out easily. Writing the Fourier component proportional to e~ luJt as $, let us consider the 
equation of the following form, 

^ 2 $ = A^> = ("J^F + F ( r *)) ( 34 ) 

where the operator A is the self-adjoint differential operator and V(r*) is a smooth function of a coordinate r*. (As 
we see later, r* corresponds to the tortoise coordinate, dr* = f~ 1 dr.) Then, if the operator A with domain C^°(r t ) 
is a positive symmetric operator in the L 2 -Hilbert space with respect to the inner product 

(Sx,^)^ = / dn §i*(r*)Mr*), (35) 

the system does not have normalizablc negative mode solutions. Consequently, the amplitude of the solution remains 
bounded for all times as long as a smooth initial data of compact support in r* is concerned (See [HI for the 

argument of initial data.) 

We should notice that this stability condition of positive self-adjointness is not a necessary condition, but is just a 
sufficient condition in general. In fact, for some type of potential which is not positive definite, it is possible to prove 
stability of the system by shifting the bottom of potential. The method is known as S -deformation: Introducing a 
new differential operator 

D=-^ + S( rie ) (36) 

with S being some function of r*, the inner product is evaluated after integration by parts as 

($,^$) i2 = Jdr* (|/3$| 2 +7|$| 2 ) , 

V = V+^-S 2 . (37) 

where the boundary term vanishes for $ £ C£°(r*). Therefore the S-deformation shifts the bottom of potential. 



B. Tensor perturbations 



The master equation JjJj for the tensor- type perturbation is the same form of Eq. I|34|) . Fourier-expanding along z 
direction, the operator A is given by 



A=-^a+fkl + V T . 



(38) 



where r* = J drf^ 1 ^) and k z is the wave number in z direction which corresponds to the mass spectrum of Kaluza- 
Klcin (KK) modes on z = z c plane. The mass spectrum is taken to be fc 2 > without loss of generality. Otherwise 
the linear perturbations break down at some z, even at an initial time, fc 2 > is called massive modes, and fc 2 = is 
zero-mode which corresponds to the perturbations of the higher dimensional Schwarzschild black holes. 
For the background given by l|33|) . the potential Vt is expressed as 



Vr(r) 



f 



n(n + 2) n(n + 1) 



( r / r h) 



n-l 



k\ - (n - 2) 



(39) 



3 It is interesting to study the stability of the black string in the other backgrounds, for example, with K ^ 1 and the cos mological 
constant. However, for such cases, even the stability of the Schwarzschild black hole has not been established completely llSl . 
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Since the spectrum of kj, satisfies fcf, — (n — 2) = (I — + n) > 0, the potential Vp is positive definite in the 
Schwarzschild wedge. Therefore we conclude that the black strings are stable with respect to tensor perturbations. 

This result is easily understandable. The operator (|38|l is nothing but the same one as the higher dimensional 
Schwarzschild black holes, except the presence of KK modes. The massive modes increase the stability of perturbations 
due to its positive contribution. This completely fits in with our physical intuition, and it might be anticipated that 
other type of perturbations are also stable due to the massive spectrum of KK modes. However, as we see below, 
the master variables of vector and scalar perturbations for the zero mode cannot be used as master variables for 
massive modes. The massive modes give new physical degrees of freedom for vector and scalar perturbations and the 
transversely gauge-invariant equations give coupled second order differential equations. Then the naive expectation 
like the tensor perturbation does not hold. 



C. Vector perturbations 

1. Stability of t = 1 

The equation for t = 1 is given by (|19|l . By introducing a new variable F W = r"™/ 2 *, we can transform the 
equation into the form of l|34|) with potential 



2(1 - n) + (2 + 3n)f . (40) 



V ~ r 2 

However, the form of potential V v is not positive definite. It becomes negative near the horizon for k 2 r 2 < (n 2 + 
n — 2)/2, and the stability for such light modes are not obvious. 

The positive definitcness of the symmetric operator A with the potential l|4l)|) is shown by the S'-dcformation. We 
find that the following choice 

gives positive definite potential V = k 2 f > 0. Therefore, this mode which corresponds to adding a rotation to the 
background solution is dynamically stable. 



2. Stability of £ > 2 



Instead of solving Eq. (|18J) . which gives coupled differential equations, let us consider Regge- Wheeler gauge by 
taking Ht = on z = z c surface. In this case, the dynamics of F a = f a arc given by which in general gives two 
coupled differential equations. We introduce the following new variables after Fourier-expanding F a in z-dircction. 



From l|lt>| l. ^ is solved as 

and we find an equation for $, 
□ - 
V v = 



Vv 
f 
f 



$ = 0, 
1 

1 + TZ 2 



F 1 



\J my + k 2 r 2 

[nm v + (n + 2)k 2 z r 2 ] $ + 2r(m v + k 2 z r 2 )& 
2(m y + k 2 r 2 ) 3 / 2 



dt + h(r), 



(42) 
(43) 



(44) 



(45) 



(n-2) + 



3/ 



1 +TZ 2 



(l+K 2 )m v + 



l-2(n + 2)/ 



{1 + K 2 ) 



(n + 2) 



/(3n + 2) - 2(n - 1) 
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where 1Z 2 = k 2 z r 2 /my. Here t„ is an initial time and h(r) is an arbitrary function. The Eq. p6|l contains only first 
time derivative of F*, and hence the initial data of F t can be specified only by h(r). 4 

The potential Vy becomes negative near the horizon for «>1. However, the positive definiteness of this potential 
can be shown by employing the S'-deformation. Applying the 5-dcformation (|41fl . the last term in the curly brackets 
are cancelled out. Using the fact that for K = 1 and I > 2, my is bounded below as my >n + 2, and then the second 
term is easily shown to be positive definite. Therefore, we conclude that the vector perturbations are stable. 



IV. SCALAR PERTURBATIONS 



A. Gregory-Laflamme mode 



The s-wave (l = 0) perturbation is the unstable mode studied by GL. Here, we discuss this mode in our framework 
and recover their result. We can use the residual gauge degrees of freedom l|A7|l to fix unphysical gauge modes. 
After eliminating the terms proportional to z-derivatives of F\ and F by using Einstein's equations, we can apply the 
harmonic gauge condition (|32|) to rewrite F\ and F on z = z c in terms of F[ and F£ . Then we finally obtains a second 
order ordinary differential equation (ODE) of F£ (or F[) in Fourier space, assuming F£ oc e ot + lfc ^. Although it is a 
second order ODE with respect to r, the equation in the original space contains higher derivatives of t and z. (See 
|36l | for more tractable equation.) From the master equation, the boundary conditions required for a normalizablc 
mode are 



F r K J^Z rfc )l-n/(n-l) ■ ( r ^ r h) 



(46) 



This is a one-parameter shooting problem with shooting parameter f2 > 0. We have solved this problem numerical! 
searching for the growth rate fl for given k z . The result is shown in Fig. Q which agrees with the original analysis 

5 



ally, 

m 



Another type of simple master equation can be obtained by taking static limit. To obtain the static mode, we adopt 
a gauge fixing 



without fixing the pure radial gauge T r (r) 



F tr = 0, (Z = Z c ) 

In Fourier space, we find a master equation 



(47) 



" d 2 


1 










dr 2 







n 



0. 



V = 2N A + N(N +7N + 12) f 6 - 2(3A^ + 6A^ - 2k 2 r 2 )f 2 - N(ZN Z + N - Ak 2 z r 2 )f 
Q = A(V 2 (2A 2 + 6N - 3k 2 z r 2 ) + N{2N 3 + 6N 2 - 3k 2 z r 2 N - 8k 2 z r 2 )f 



-2/(2A 4 + 6A^ 3 



Zklr 2 N 2 



Ak 2 r 2 N- 



2 k y) 



K = r 2 f [N 2 + N(N + 4)/ 2 - 2(iV 2 + 2N - 2k 2 z r 2 )f] 



(48) 



where N = n — 1. Other components F£ and F are given in terms of F\ . (Note that another type of master equation 
was derived in |37|. which is more tractable than the above equation in practice.) A Neumann condition on the 
horizon is obtained by requiring the regularity on the horizon. Solving this equation is the one-parameter shooting 
problem with the shooting parameter k z . Hence we can think of this equation as a master equation determining the 
GL static mode. As is known well, the wave number of this static mode, which will be denoted fc C rit, gives a critical 
point at which stability of the s-wave perturbation changes. For k z < k CI it, the perturbations are unstable, whereas 
they becomes stable for k z > fc cr i t (Fig. ^ . 



4 There will be another arbitrary function. Substituting the solution of 1451 into the Eq. 1151 . we can integrate it by z to get Hx at 
z z c . Two arbitrary functions of x a appears, but one of them can be eliminated by The remaining function corresponds to an 
"initial data" in the bulk, whose evolution is stable sine it is homogeneous (zero mode) in ^-direction 

5 The harmonic gauge does not fix the gauge T a completely. Besides the static radial gauge transformation, the residual gauge is Tt oc r 1— " . 
F[ depends on this gauge while F£ is free from this mode. 
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FIG. 1: Plot of fl as a function of fc z for black strings with spacetime dimensions D — n + 3 = 5, 6, ■ ■ • , 14. The wave number 
k z is normalized by the horizon radius rv 
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FIG. 2: Static mode search. The possible asymptotic solutions of Ft are Ft cx e ±kzT . In the figure, we plot Ft at some 
r/rh 3> 1 with respect to the single shooting parameter k z . Since the normalizable solution decays exponentially, each narrow 
"throat" corresponds to a normalizable mode. The critical wave numbers agree precisely with the static limit f2 = in Fig. Q 



The equation l|48|l has two asymptotic solutions behaving F\ oc e ±fcr , and only the decaying mode is the physical 
normalizable solution. Such physical solution can be easily found by searching a minimum value of F\ at some fixed 
asymptotic point as a function of k z . Figure |21 shows the result of the shooting problem. As we see, the critical wave 
number fe cr it can be precisely determined by this method, and this agrees completely with the analysis of dynamical 
perturbations discussed above. 



B. 



For the zero mode k\ = 0, the 1=1 mode has no physical degrees of freedom. This can be easily observed from the 
fact that the master variable of the massless mode can be reintroduced by recovering the lacked equation l|C2|) as a 
gauge condition. However, since such gauge fixing is not complete, there remain additional residual gauge degrees of 
freedom. By using the residual gauge degrees of freedom, it is shown that there is no dynamical degrees of freedom 
in the vacuum case ^| 0] . (More direct counting of physical degrees of freedom is possible by taking F = gauge 
fixing.) 

For the KK modes, there is no unstable dynamical mode. This is confirmed directory by solving the EOMs on 
z = z c . Let us take the gauge F = 0. This is not a complete gauge fixing, but F a b does not depend on the residual 
gauge. After eliminating all terms proportional to d 2 z F by employing (|30a|l . we can solve the EOMs explicitly after 
tedious calculations. One finds that only trivial solutions are allowed on z = z c in the present case, so that there is 
no unstable dynamical degree for the KK modes. 6 



Equation I30cl for the KK modes with 1311 corresponds to the constraint equation ICll for the zero mode. Eq. !3Ual is used to take 
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I > 2 



For the generic modes of scalar perturbations, the Einstein equations consist of five equations, and they give coupled 
partial differential equations. We decompose F and F% as follows; 

(V> - $P - q) 

1 . 

^{tp-p-5q) 
(ip + 3p + 3q) 



F = 


i 


gr n-2 


Ft = 


1 




f; = 


i 






-d t Z 

r n — 2 



(49) 



Then from the transversely gauge- invariant equations (|24[1 and (|25|) , we obtain 

„ 9 kl (l-n+(n + S)f) 2(1 -n + nf) (nf-2n + 2) , 
□Z + d\Z - -f Z + V JJ p + — ^—^q + — — Z' = 0, 

iyt£t ij* IJ* ^ 



(50) 



k s + 



4 - 4n + (3n - 10)/ 



□ g + ^ + 2 ^-^ 1 + 3 ^ z-f 

3r/^ 3r 



2(1 - n) + (3n - 2)/ 



3r 2 3r 



(3n - 4)p' + iq' 



= 0, 



8p' + (3n - 8)q' 



(3n 2 + 6n - 25)/ 2(n - 1) (n - l) 2 



2/ 



(3?i 2 - 9n + 5)/ 2(n - 1) (n - l) 2 



/ 



— fc| 



and ip is given by 



( n + !) «a. 



(1 + 5n) A2 „ , (5 ^ _ _ x + + p; _ _ x + ✓ 



p_ 

2^ 

2 



(n 2 +3n-4)-(n-5)/- 



2p + q 



[n-lf 
" ;d?d r 



2q 



y/fZ 



(n 2 -3n + 2) + (2n-l)/ 



(n-1) 2 



f3/2 t 

From other remaining equations, we obtain a non-trivial equation for ip 
2 [1 - n + (n - 3)/] - , V 

(n-1) 2 



Z + £ [(„ - 2ff + (n - 2)(1 - n) - *§] + - ^(2p' + <?') 



(9n 2 - 56n + 71)/ - 2(n - l)(5n - 16) + 



/ 



2r 2 



(lln-17)/ + (n-8)(n-l) 



(n- l) 2 



/ 



(51) 







We first notice that Eq. (j50|) does not contain tp, and i/j can be determined by (|51|) once we solve p, q and Z. Hence 
it is sufficient to analyze Eq. H50|) for the stability problem. We begin with a limited case to study the stability. In 
the limit ks > 1, the EOMs are 



□ 



□ - \ki 



k 2 



■/■- 



0. 



p = — z, 

3rf 



kg 



n 
r 



q = 



{n-lf 2[3(n- 1) - (3n + l)/] 



2r 2 / 



P 



3r/ 2 



(52) 
(53) 
(54) 



where we have left the terms proportional to 1// since it becomes dominant near the horizon. Apparently, Z are 
stable due to the positive definite potential. Then, the stability of p and q is also obvious. Furthermore, we notice 



the gauge F = 0, and Eq. I30bl works as a constraint equation. 
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FIG. 3: Stable-unstable phase on k\-k\ plane. The zero mode k 2 = corresponds to the perturbations of higher dimensional 
Schwarzschild BHs, which are stable. The Gregory-Laflamme mode (I = 0) is at k% — with k 2 < k 2 lit . On the plane, the 
shaded upper-right corner with ks S> 1 or k 2 3> r^ 2 is shown to be stable analytically. The stability of other generic modes is 
confirmed numerically. 



that the same argument holds true for very massive modes k 2 3> r h 2 , without taking the limit » !■ Thus the 
system is stable k$ ^ 1 or k 2 ^> r^ 2 . We note that the system is stable in the zero mode limit k 2 — > 0, since in 
this limit the perturbations are the same as the Schwarzschild black holes. Hence on the k 2 - fc| plane, there exists 
stable region. The stable/unstable parameter region discussed here is summarized in Figure [3] 

For general modes with arbitrary k 2 and kg, we performed a numerical search for unstable solution, as we do 
in Sec. IIV A1 Assuming an unstable perturbation oc e nt e lk ^ z 7 we obtain boundary conditions similar to (|46|l . Then 
we performed a parameter search in the relevant region of (fc z ,f2) and no solutions were found, suggesting that no 
instability exists for the generic modes. To confirm this result furthermore, we have also performed a search for 
critical static mode: if the system is unstable, a static mode will exist since the real eigenvalue in the stable region 
will cross the zero axis at lease once when it becomes unstable. Since the horizon boundary conditions are not the 
same as those obtained by just taking the static limit of dynamical perturbations, this numerical search works as an 
independent search of unstable mode. Redefining dtZ = we take static limit of <|5(J[) . In this limit {p, q} and £ are 
decoupled, and we can easily performed the search. The differential equations for {p, q} are a two-parameter shooting 
problem, and a part of the result is shown in Fig. 0] which corresponds to Fig. [3 Clearly, there is no static solution 
satisfying appropriate boundary conditions. The same result holds also for £. Therefore we conclude that the black 
strings are stable for all types of perturbations except the s-wave mode. 



V. SUMMARY AND DISCUSSION 



In summary, we have studied stability of black strings with respect to all types of gravitational perturbations. There 
are three types of perturbations; tensor, vector, and scalar perturbations. The vector and scalar perturbations have 
the exceptional modes of multipolc moment besides the generic modes. For the higher dimensional Schwarzschild 
black holes, the exceptional modes are not dynamical degrees of freedom. However, we have paid particular attention 
to the exceptional modes since they might become dynamical with some instability. 

The generic modes of tensor (£ > 1) and vector [I > 2) perturbations have been shown to be stable. The generic 
modes of scalar (£ > 2) perturbations were studied partially employing numerical investigation, and they have been 
shown to be stable. For the exceptional modes, we have discussed that the vector perturbation of £ = 1, which 
corresponds to adding a rotation, is stable, and the exceptional mode £ = 1 of scalar perturbation has no unstable 
dynamical degree of freedom. The £ = mode of scalar perturbation is also the exceptional mode, and it is dynamically 
unstable as discussed by Gregory and Laflammc. After all, the unstable mode of gravitational perturbations for black 
strings is only the £ = mode of scalar perturbation. 

The zero mode (k\ — 0) of the scalar perturbation with £ = corresponds to a shift of the mass parameter of 
the higher dimensional Schwarzschild black holes (or uniform black strings), and hence this mode is not allowed as a 
consequence of the Birkhoff 's theorem. However, the KK mode with £ = is essentially different from the gravitational 
perturbations of the Schwarzschild black holes, and in fact it does not change the mass of the black strings. Therefore, 
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(dp/dr) 



FIG. 4: Search for critical static mode for I = 2, 3. (See Fig. [21 for I = 1 mode.) This is a two-parameter shooting problem. 
The two parameters are the wave number k z and the derivative of p at the horizon. The figure shows a plot of (p 2 + q 2 ) at 
some asymptotic region with respect to the two parameters. Possible asymptotic solutions are p, q oc e ±fc ~ r , and normalizable 
solutions will decay at r > rv No narrow "throat" appears so that there is no normalizable static mode. For other higher 
multipoles (£ > 2), we obtained the same results. 



from the viewpoint of effective theory on a z = const, plane, we understand that the existence of Gregory- Laflamme 
instability is directly related to the inapplicability of Birkhoff 's theorem. 

This observation is useful to consider a possible counterexample of correlated-stability conjecture (CSC). If we do 
not interpret CSC in strong sense, the instability predicted by CSC is s-wavc instability |38j . To have some insight, 
let us discuss a black hole obtained by dimensional reduction of the black string/brane. If the black hole is a hairy 
black hole and (generalized) Birkhoff's theorem cannot be applied, the s-wave perturbation becomes a dynamical 
degree of freedom. This s-wavc perturbation is homogeneous (zero-mode) perturbation in the original spacetimes, 
and it is not the (massive) perturbation for which CSC concerns. Then if there is a model in which the homogeneous 
s-wave becomes unstable for some parameter region, the model will be a counterexample of CSC since the instability 
is disconnected from CSC. In fact, the recently proposed counterexample is based on a hairy black hole H^,^3' an d 
the unstable mechanism is along the line of the above discussion. 

We finally address possible extension of present analysis. First, it is interesting problem to study the stability of 
charged black strings, focusing on how a given charge works to make the string stable near the thermodynamically 
stable and/or BPS state. It will give us deeper understanding of CSC from the perspective of dynamics. Second, 
we have analyzed the black object with a single trivial transverse direction, for simplicity. For black branes with 
translationally invariant multiple directions, it will be possible to expand the perturbation variables by harmonic 
tensors associated with the uniform transverse directions. We would like to discuss these issues somewhere else. 



Acknowledgments 

It is a pleasure to thank Kcnta Kiuchi for helpful correspondence and discussion at an early stage of this project. 
The author would like to thank Shinji Mukohyama, Hideo Kodama, Akihiro Ishibashi, Umpei Miyamoto and Robert 
C. Myers for their very helpful comments, discussion and conversation. In particular, Akihiro Ishibashi presented a 
useful observation. He would like to thank the organizers of the "Scanning New Horizons: GR Beyond 4 Dimensions" 
workshop in Santa Barbara for a stimulating environment, a relaxed atmosphere and hospitality while this work was 
being completed. The author also wants to thank Barak Kol, Dan Gorbonos and Hebrew University of Jerusalem for 
their hospitality and stimulating discussion during an initial stage of this work. This work was supported in part by 
JSPS (Japan Society for the Promotion of Science) fellowship and the National Science Foundation under Grant No. 
PHY99-0794 



14 



APPENDIX A: GAUGE TRANSFORMATION 

In this appendix we summarize the transverse gauge transformation Ipjl. The metric perturbation Hab transform 

as 

Sh AB = -VaZb - VbU (Al) 

in terms of the infinitesimal gauge transformation dx A = £, A . The transverse gauge transformation JSJ) can be 
decomposed into 

Sh ab = -Dah-Dbta (A2) 

6 



5h ai = -r z D a l^j - Diia (A3) 
Shu = -D i ^-D j Ci-2r lij ^ a D a r. (A4) 

Since the infinitesimal transformation £ has no tensor component, the expansion coefficient of the tensor perturbation 
is gauge invariant. Our interest is therefore gauge transformation of vector and scalar perturbations. 
The vector component of the transverse gauge transformation is 

£ a = 0, 6 = rLVi (A5) 

for the modes k v ^ (n—l)K, where L = L(x a ) is an arbitrary function. Then the corresponding expansion coefficients 
of the perturbation transform as 

6f a = -rD a (-), SH T = ^L. (A6) 
\r J r 

As for the scalar perturbation, the gauge transformation for k 2 s (k 2 s — nK) ^ are given by 

£ a = T S, fc = rLSi. (A7) 

Under these transformations, the expansion coefficients of the metric perturbation transform as 

5 fab = -D a T - D T a , 
5 f a = -rD a ( -) + ^T a , 



r J r 

ah ks r^ Dar T 
oH L = L T a , 

nr r 

SH T = —L. (A8) 
r 

The gauge transformation for kg(kg — nK) = are obtained by setting appropriate functions equal to zero in the 
above equations. 

APPENDIX B: DETAILS OF CALCULATIONS 

In this Appendix, we summarize the details of calculating perturbed Einstein's equations for completeness. Some 
of them are based on Ref. [U ■ 

1. Background Quantities 

We consider perturbations of spacetime on (n + 2 + l)-dimensional spacetimc whose unperturbed background 
geometry is given by the metric Jp). Decomposition of connection coefficients is 

n c = {2) n c (y), r?- = -r(D"r) 7 «, = f U = r^O*)- ( B1 ) 
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Here ^ 2 ^T bc is the Christoffel symbol of the two-dimensional orbit spacetime. Curvature and Ricci tensors are 



R" 



bed 



{2) R a bcdl 



ajb 



D a D b r 



Rat = {2) Rab - -D a D b r, & 



□r 



R i jkl = [K-(Dr) 2 ](gi lj i-gh jk ). 
K -{Dr) 2 ' 



+ (n - 1) 



■9/ 



i? = (2) i? - 2ti — + n{n - 1) ^ — '—. 



(B2) 



Einstein tensors are decomposed as 

G a b = {2) G ab - -D a D b r - 



n(n - 1) K - {Dr) 2 n 



9ab 



_ 1 (2) _ {n-l){n-2)K-{Dr) 2 n-1 
2 2 r 2 r 



□r 



G M = 0. 
For the two-dimensional metric 



(B3) 



ds 2 



-f{r)dt 2 



fir) 



-dr 1 



Ricci tensor and Ricmann tensor are explicitly given by 



(2)d til ub _ xb r> I s^R 

U — — / > -Kq — a~ 2~ ' -ftabed — \9ac9bd - 9ad9bc)—^- , 



and non-vanishing Christoffel symbols are 



2f 



■pi J "pr J J "pr 



2 ' 



2/' 



(B4) 



(B5) 



(B6) 



2. Perturbations of the Ricci Tensors 



We consider metric perturbations under the gauge fixing of Eq. ©. In general the perturbation of the Ricci tensor 
is expressed in terms of Hmn — &9mn as 



26R MN = -V L V L /l M A'-VMVAr/7 + VMV L ^ + VArV L /l^ 



5R 



-RMhh N + RMhh M — 2RMLNsh , 
-h MN R MN + V M V N h MN - V M V M h. 



(B7) 



Here and hereafter the trace h is h\ = Hmn9 = + r h l3 lij 



a. Decomposition formula 



To calculate the perturbed Ricci tensor, we need to decompose the connection V into D and D. The operator D 
and D work as 



Djh ab 

Djhai 
Dahtj 

D a h bj 



= djh ab , 

= djhai - 

= d a h bj -Wr e ab h ej . 



1 ji'l-ak, 



(B8) 
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The followings are useful formulas of decomposing the operator D for arbitrary tensor h A B and vector T A . 

W a T b = D a T b , 



and 



Vjlj = DiTj + r{D a r)^ijT a , 

ViT a = DiT a - — Ti, 
r 

V a T, = D a T l - —T,, 
r 

V z Ta = d z T A (B9) 

V a h c d = D a h c d, 
V Q Vfcft- C( j = D a Dbh c d, 
D a r 

^ ahij D a hj i j 2 H^j 3 

V a hbj = D a hbj ~hbj, 

D b r D c r 

Vi/ibc = Dihb c hi c hbi, 

r r 

D r 

Vihjc = Dihjc + r(D a r)^ijh ac —hj, 

Vih, jk = Dihjk + r(D c r)jijhck + r(D c r)ji k h JC , 
ViVjh = DiDjh + r(D c r)j i:j D c h, 

V z h AB = d z h AB . (BIO) 

b. Perturbed Rica tensor 



26R ab = -Uh ab + D a D c h% + D b D c h c a 
D c r 



-n (-D c h ab + D a h cb + D b h c 

r 

- {2) R a h cb + i2) R b hca - 2 (2) Racbdh 



cd 



\{D a D l h bl + D b D l h ai ) - -^-D a hi^ 



;Ah ab 



—D a rD b rh in l i - D a D b h - d z z h 



a b • 



(Bll) 



25R ai = D t D b h a 



b , n —±D h rD i h, 



i'l>ab 



-rD [ -h a 

r 



-D»rD b h ai - D a rD b 



71+ 1 



D b rD a h bl +rD a D b { -h b 



+ 




f (n-1) — 


- (Dr) 2 


Or' 










r 




+ -^D b rD a rh bt + 


(n + l)rD a 


(>) 


hbi 




n + 2 D a D b rh lb 
r 


+ {2) R b a hu- 




+ ^Dibiflaj 



rD n 



—D J hji 



-^D a rb J hji - ^D a rD l h jk j 
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-rD a [ -Dih ) - d 2 z h au 



(B12) 



25R t 



[2rD a rD b h b a + 2(n - l)D a rD b rh ab + 2rD a D b rh ab ] 7ij 



rD t D a [ -h a 3 



rDjDa 



-K 



D a r - - D a r - 

+(n - 1) — (A/iaj + £>A*) + 2 — D kh kali 



-~2 A/li 



2(j kl h k aij - hij) 



- n—Dahij + \(DiD k h kj + DjD k h ki ) 
(n-l)- + 2 l 



r r 
A" - (£>r) 2 



,2 



,(£>r 



-ji 3 ~f kl hki 



-DiD 3 h - rD a rD a h ll3 - d 2 z h tj , 



(B13) 



2n 

&R= D a D b h ab + —D a rD b h ab 
r 



_(2) R ab + ^ D a D b r + n(n ^) D a rD b r j ^ 



D a r 



-—D a D l h1 + 2{n - l)^-D l h 



1 - ,■ - , D a r • ■ 1 

D a r 1 ~ 

-□/i - n D a /i - — A/i - <9|/i, 



A (£>r) 5 



( n _l) 2 



SRaz components are 



28R lz = d z {-£>ih + D b h\ + ^b 3 h 3l H 



D a r 



r 

D c r 



— hci\ 
r ) 



25R Z 



-d 2 h 



(B14) 

(B15) 

(B16) 
(B17) 



APPENDIX C: EINSTEIN EQUATIONS 



Einstein equations for scalar perturbations are summarized as follows. From the components SGf and traceless 
part of SGf of the Einstein equations, we find the following equations: 



k s 
2r* 



L-D b (r n ~ 2 F b )~rD/ 



2{n-2)F + F' c 



2(n-l)D a F 



+ rf a .. 



0, 



(CI) 
(C2) 



6G ab and 5G\ gives another two equations. 

- hdG ab = DF ab - D a D c F bc - D b D c F ac + n—(D c F ab - D a F bc - D b F ac ) - R^ c a F cb - R^ c b F c , 



+ 2 RaL FCd ~ ^§ p ab + D a D b F c c + 2n \ D a D b F + X -D a rB b F + X -D b rD a F 
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D c D d F cd + —D c rD d F cd + ( —D c D d r + n(n - X) ElI^L _ fl (2)cd ) Fcd _ 



2n(n+l) 



D c rD c F + 2(n - 1)^— ^F - DF C C - -D d rD d F c c + %F r c - / C C . Z2 - 2nff£, 



(C3) 



1 1 

StT 



1 fr; _ 1 1 n a r 
> b F b - -OF: + { - J- (2D b F b a D a Ft 



2r 

-(n-ipF-^fc^^W + %^ 1 
r 2nr z 

Explicit equations from SR z a — components are 

25R ZZ = -Sd 2 z (fc + 2nH L ) = 0, 
2d z 5R az = S 



2r 2 



2{n-2){k s -nK)F + k s F« - - (n - 1)H L , ZZ (C4) 



(C5) 
(C6) 



2&£R« - 8, flg £» c (r/ c ) + n(£) c r)/ c - 2k s H L + 2H T 



n — 1 kg — nK 



n ks 



0, 



(C7) 



where we have used <|C5(1 in (|C6|) and l|C7(l . 

Let us try to rewrite Eqs. (|C3|) and l|C4(l . Taking the trace of (|C3|I and combining it with (|C4(1 . we can solve 
and Hl, zz in terms of _F and i 7 ^;,: 



2(n + l)F L , 



UF C C - D a D b F b - —(D a F b b - 2D b F b ) 

kl ™ 2n(n + l) fl j, , 2(n-l)(n + 2)*|-nK' 



1 .Da.Dbr 



2(n + l) + (n- l)(n + 2) 



F 



D a rD b r m 



(n + l)f a a , zz = -2nUF c c + 2nD a D b F ab - n(n - \)^-{D a F b - 2D b F b ) + (n - 1)^J F c c 



D a rD b r 



F ab - 2n(n - 1)DF - 4(n - 1) 



ki-nK 



-F 



ab 



(C8) 
(C9) 



Substituting these into 1|C3(I we obtain 



UF ab - D a D c F bc - D b D c F ac + n —(D c F ab - D a F bc - D b F ac ) - R%> F b c - R { $ F* 
+ 2 Racbd FCd - % p ab + D a D b F c c + 2n ( D a D b F + -D a rD b F + -D b rD a F 



9ab 
71 + 1 



D c D^ cd - ^D c r(D c F d - 2£) d F c d ) + + n(n 1)^^) F cd 



2n(n + l) r , ,kl-nK 
-2nUF v D rD c F + 2{n - 1)^—5 



F - DF5 + ^-F, c 



r 2 " c 



+ fab,zz = 



So far. we have only used (|C3|I and l|C4(l . From IjClf) and (|C2() . one can construct X a as 
Then using this, we can rewritten (|C10(I and (|C8|) in terms of F ab and F, resulting in Eqs. (|2"4*j) and l|2"5)l. 
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